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We investigate topology in lattice simulations of QCD with 
two flavours of dynamical Wilson fermions. At various sea 
quark masses we find reasonable agreement between results 
for the topological charge from fermionic and gluonic defi- 
nitions, the latter with cooling. We study correlations of 
the topological charge with light hadronic observables and 
find effects in the flavour singlet pseudoscalar channel. The 
quark mass dependence of the topological susceptibility, Xi 
is consistent with the leading order small m, expectation, 
X = 

PACS numbers: 11. 15. Ha, 12.38.Gc, 11.30.Rd, 12.38.Aw 



I. INTRODUCTION 

One of the most intriguing features of QCD is its topo- 
logical vacuum structure which results in phenomenolog- 
ical consequences with an important bearing on particle 
physics such as the breaking of the axial Ua(1) and chiral 
symmetries. 

In the context of pure gauge theories, lattice meth- 
ods have by now matured to provide a valuable tool 
kit for accessing the continuum topological susceptibil- 
ity, x- Various gluonic definitions of the topological 
charge lead to consistent estimates of x |), and rem- 
nants of the Atiyah-Singer index theorem are witnessed 
to hold on the lattice: the topological charge, Q, as 
obtained after cooling from the gluonic degrees of free- 
dom is consistent with the fermionic index, as determined 
from counting the zero-level crossings of the eigenvalues 
of the Wilson-Dirac operator, 75(^4- too), under varia- 
tions of the bare quark mass, too ~ MM. Moreover, 
one might perceive the very value of x as delivered by 
quenched lattice simulations x\ /A = 213(14) MeV, 

to confirm the lattice approach to topological proper- 
ties: it turns out to be in fairly good agreement with the 
large N c anticipation made by Witten and Veneziano |S , 
X q w f%(m*, + m\ - 2m 2 K )/6 « (180 MeV) 4 . In conclu- 
sion, the lattice machinery appears to work for the study 
of topological aspects in gluodynamics. 

Unfortunately, however, the situation regarding the 



QCD vacuum proper is by far less settled: for sheer cost 
reasons so far only the regime of intermediate sea quark 
masses, not much lighter than the strange quark, has 
been explored. In fact it is still under debate whether 
the data from full QCD simulations support the ex- 
pected vanishing susceptibility at chiral sea quark masses, 
to — > 0: 



X = — S + C(to 2 



lim lim (ipif>). 

m^0 V^oo 



(1) 



Attempts to verify this prediction via lattice simulations 
have a long history. While the first lattice studies in this 
direction with Kogut-Susskind (KS) sea quarks |-|| 
could not reach conclusive results, later investigations did 
yield some qualitative evidence in favour of the expected 
decrease of x( m ) @- 

Quite recently, this problem has been revisited from 
three different sides — but the debate is still open: the 
CP-PACS U and UKQCD || collaborations employed 
improved Wilson fermionic and different gluonic actions 
while the Pisa group ]l6| operated with two and four 
flavours of KS fermions. A further data point for nf — 2 
KS flavours has been obtained by Hasenfratz JlTj (for 
a recent review see e.g. Ref. |lql). While both the Pisa 
group (working at an inverse lattice spacing a -1 ~ 2 GeV 
and bare quark masses otoo > 0.01) and the CP-PACS 
collaboration (at a -1 w 1.3 GeV, m^/mp > 0.59) saw 
no evidence whatsoever in favour of the expected chi- 
ral behaviour, the UKQCD collaboration did verify (at 
a -1 ss 2 GeV and m^/nip > 0.57) a decrease of x( m )> 
consistent with theoretical expectations. 

In this paper we shall present an analysis focused on 
this issue, based on the final statistics of our SESAM and 
TxL samples of QCD vacuum configurations. The arti- 
cle is organised as follows. In Sect, [n] we describe the 
details of our simulation and methodology, including a 
comparison between fermionic and gluonic definitions of 
the topological charge. Preliminary results on this com- 
parison, based on smaller statistical samples, have been 
reported by us previously fll9| |. In Sect. [II we investi- 
gate correlations between the topological charge and the 
hadron spectrum. Finally, in Sect. IV, we present the 



lattice data on the topological susceptibility. 
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TABLE I. Simulation parameters. Unless an error is stated 
the autocorrelation times T int are lower estimates. The num- 
ber of measurements is related to the number of trajectories 
by n moaa = iVeaa^tiaj , with ^ mcaa = 1/10 in the quenched case 
(J3 = 6.0), i/ me a 8 = 1/24 at ft = 0.1565, v mMS = 1/2 in the 
second series of measurements at n — 0.1575, L a = 16 and 
= 1/25 otherwise. The physical units have been ob- 
tained by setting r^ 1 — 394 MeV and are subject to scale 
uncertainties of about 5 %. All unquenched simulations have 
been performed at /3 — 5.6. 



K 






rod 1 


a'i/GeY 




Tint 


0.1560 


16 


206 


5.11(3) 


2.01(1) 


7.14(4) 0.834(3) 


30 


0.1565 


16 


209 


5.28(5) 


2.08(2) 


6.39(6) 0.813(9) 


55 


0.1570 


16 


200 


5.48(7) 


2.16(3) 


5.51(4) 0.763(6) 


40 


0.1575 


16 


215 


5.96(8) 


2.35(3) 


4.50(5) 0.692(10) 


55 






2200 


53 


11 


11 11 


54(4) 


0.1575 


24 


150 


5.89(3) 


2.32(1) 


6.65(6) 0.704(5) 


140 


0.1580 


24 


140 


6.23(6) 


2.45(2) 


4.77(7) 0.574(13) 


130 








6.73t[i 


2.65+| 


0.179(1) 




P = 6.0 


16 


320 


5.33(3) 


2.10(1) 
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II. DETERMINATION OF THE TOPOLOGICAL 
CHARGE 

A. Measurements 

We analyse ensembles of gauge configurations that 
have been generated by means of the hybrid Monte Carlo 
(HMC) algorithm using the Wilson fermionic and gluonic 
actions with n/ = 2 mass degenerate quark flavours at 
the inverse lattice coupling, /? = 5.6, corresponding to 
an inverse lattice spacing a -1 = (2.65^g ± 0.14) GeV at 
physical sea quark masses. This was done on L 3 x L T = 

16 3 x 32 as well as on 24 3 x 40 lattices at five different val- 
ues of the sea quark mass parameter. The corresponding 
chiralities can be quantified in terms of m^/m^-ratios, 
ranging between 0.834(3) and 0.574(13) @. The rel- 
evant simulation settings are displayed in Table Q pT[ . 
At each K-value 4,000-5,000 thermalized HMC trajecto- 
ries have been generated. In addition to the dynamical 
quark simulations, quenched reference measurements on 

16 4 lattices at (3 = 6.0 were performed. The configu- 
rations that enter the analysis are separated by 25 (24 
at k — 0.1565 and 10 in the quenched simulation) HMC 
trajectories. 

Lower limits on the integrated autocorrelation times, 
Ti n t, arc estimated by binning the data for the topolog- 
ical charge Qi , i = 1 , • • • , n mcas , into ribm blocks that 
contain m successive measurements each. On each such 
block of length m an average Qj , j = 1, . . . , nbin, is cal- 
culated: Q" 1 = 0-/m)YliLiQm(j-i)+i- We determine 
fluctuations between these bins, 



1 



nbin(nbh 



"bin 

' 3=1 



(Q?~(Q)y 



(2) 



where the average (Q) ~ is calculated on the first nt,in x 
m configurations. We then estimate the autocorrelation 
times, 



1 f^Q 2 m 

Tint « t, max ~\~Fa 



(3) 



that are included in the table. 



denotes the measure- 



ment frequency. On the 16 3 x 32 volume at k = 0.1575, 
in addition to the v mcas = 1/25 time series, we also de- 
termined the topological charge with increased frequency 
i/ mcas = 1/2, with a reduced number of 10 (as opposed 
to 60) cooling sweeps. The larger frequency enabled us 
to compute the autocorrelation time from the autocor- 
relation function itself, with full control over statistical 
errors [^9 22 1. The resulting value, T; n t = 54(4), is con- 
sistent with our estimate obtained in the way described 
above, T; n t ~ 55. We take this as an indication that our 
estimates are reasonable. Interestingly, the autocorrela- 
tion times grow both with decreasing quark mass and 
increasing volume as already observed in Ref. |l9| . 

We employ a gluonic as well as a fermionic definition 
of the topological charge. In the gluonic case we define a 
topological charge density, 



q(x) 



1 



1 



; F(x)*F(x), (4) 



Km 2 ^ v y ^ v ; 16^ 2 
where we use the symmetric (clover leaf) definition 



Fp V = ^-j (IV + U,,,-,, + U- Mt - V + U- U , M - 4) , (5) 

and F^ = he^ upa F pcr , which is correct up to order a 2 
lattice artefacts. The topological charge, 



16tt 2 



(F,*F), 



(6) 



should then approach integer values as a — > in the 
infinite volume limit, on sufficiently smooth gauge con- 
figurations. 



B. Renormalization 

In a quantum field theory, both Q and \ undergo mul- 
tiplicative renormalization. Moreover, x requires addi- 
tive renormalization. Cooling (|JJ2^] is meant to filter 
for (semi)-classical features; indeed, after cooling, both 
renormalization constants have been shown to be close 
to their trivial values [QJ|] . In the context of this article 
we will not investigate properties of the density distri- 
bution q(x) itself but rather stick to the net topological 
charge Q only. For this purpose the iterative application 
of cooling by simple minimisation of the Wilson plaquette 
action is appropriate. To keep the cooling update local 
we visit the lattice sites in an even-odd pattern, rather 
than in the sequential ordering that is usually employed 
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TABLE II. Estimates of the fermionic renormalization con- 
stant Zp, Eq. (0), assuming the gluonic one to be unity. 



K 


L, 


mo<J 


Zpuioa 


z P 


0.1560 


16 


0.0498 


0.060 


1.20 


0.1565 


16 


0.0395 


0.056 


1.43 


0.1570 


16 


0.0293 


0.045 


1.55 


0.1575 


16 


0.0192 


0.029 


1.49 


0.1575 


24 


0.0192 


0.031 


1.64 


0.1580 


24 


0.0092 


0.017 


1.79 



in the Monte Carlo updating of pure gauge configura- 
tions. The inner-most loop within a cooling sweep runs 
across the directions /i, and we carry out 60 such sweeps. 
In Fig. [I] we illustrate the resulting numerical pattern of 
topological charges which cluster nicely around integer 
values as anticipated. 
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FIG. 1. The topological charge Q determined by the glu- 
onic definition, Eq. (|^), after cooling on the 16 3 x 32 lattices 
at K = 0.1575. The results are concentrated around integer 
values (dashed horizontal lines). 

The fermionic method consists of determining the 
topological charge a la Smit and Vink [Eq], 



-Zpm tT^ 5 M 



(7) 



that is inspired by the continuum axial divergence re- 
lation, 9 M [^(a;)7 A1 75?/'(x)] = 2rmp(x)"f 5 ip(x) - 2n/g(x). 
The above trace is to be taken over Dirac, colour 
and space-time indices only (not over the flavours), i.e. 
"075^ — riftrj^M^ 1 for nj mass degenerate quark 
flavours. Zp = 0(1) is a renormalization constant, 



m 



n c )/(2a) denotes the bare quark mass 



m , 



and M is the lattice discretised version of Jf)- 
2nM xy = S xy - «2J[(1 - 7m) U x ^S x+ p, tV 

+ (l + 7/.)^- A ,/x-A, tf ]. (8) 



We determine tr75M _1 using Z2 noisy sources with di- 
agonal improvement as detailed in Ref. 27 1. On the 
k = 0.158 configurations 100 such estimates were per- 
formed while at all other k values we averaged over 400 
estimates]]. 




-6 -4 -2 2 4 6 
(F,*F)/(16ji) 

FIG. 2. Correlation between the values of Q as computed 
from the gluonic [horizontal axis, Eq. (^)] and fermionic [ver- 
tical axis, Eq. (0), Z P = 1.49] definitions on 16 3 x 32 lattices 
at k = 0.1575. The width of the error band corresponds to 
typical statistical uncertainties, due to the stochastic estima- 
tion of the fermionic values. 62 % of the data points lie within 
this la region. 

The renormalization constant Zp in Eq. (0) is un- 
known. We attempt to estimate the combination Zpm^a 
from the ratio of the gluonic [Eq. (^[)] and of the fermionic 
[Eq. (Q))] definitions. We determine m^a using the crit- 
ical value, k c = 0.15849(2) [|FJ. The results are dis- 
played in Table |l| assuming the multiplicative renor- 
malization of the gluonic definition after cooling to be 
unity. The estimates of the fermionic traces are subject 
to statistical uncertainties, AQ/Zp a 1 on the 16 3 x 32 
lattices and AQ/Z P sa 2 - 3 on the 24 3 x 40 volumes. 
In addition, both definitions are expected to suffer from 
different 0(a) lattice artefacts. We do not attempt to 
estimate the resulting statistical and systematic uncer- 
tainties on Zp. We find Zp to be of order one and to 
depend monotonously on the quark mass, with the devi- 
ation from unity increasing when approaching the chiral 
limit. 

The correlation between the two definitions is visu- 



1 Approximating Q by a finite number of noise vectors can 
result in an underestimated integrated autocorrelation time. 
Nonetheless, we find the r; n t-estimates from the fermionic def- 
inition |^2j to be consistent with the gluonic ones of Table [l| 
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alised in the scatter plot, Fig. ||, for the 16 3 x 32 lattices at 
k = 0.1575. While the gluonic data points cluster around 
integer values, this is not the case for the (stochastically 
estimated) fermionic values. The data of the figure are 
normalised such that the points should collapse onto the 
line with slope one. Indeed, 137 out of 199 values (62 %) 
lie within the one a error band. The correlation between 
the two definitions is also visible from the Monte Carlo 
histories depicted in Fig. ||. 



deed, the topological charge distribution on the large lat- 
tice at k = 0.1575 is by a factor of about two wider than 
that on the small lattice. We also observe reduced fluc- 
tuations as we increase k at fixed (3. While the 16 3 x 32 
time histories appear to tunnel ergodically through all 
topological sectors the total number of tunnellings ob- 
served for k = 0.1580 is not yet sufficient to achieve a 
symmetric distribution. 




1000 2000 3000 4000 5000 6000 
trajectory 

FIG. 3. Time histories of the topological charge Q for 
the gluonic [(F,*F), Eq. (^)] and the fermionic [tr7sM _1 , 
Eq. (pj)] definitions. The last two series have been obtained 
on 24^ x 40 lattices, the first four on 16 3 x 32 volumes. 



C. Ergodicity 

In Fig. [| we survey the Monte Carlo histories of Q 
for all our runs. These charts provide a gross evidence 
for the quality of our data in view of the decorrelation 
with respect to the topological sectors. The topological 
susceptibility, 



= £<g(x)«(0)> = 



V 



(9) 



should be independent of the volume V — L\ x L T a 4 to a 
first approximation. Therefore, the modulus of the topo- 
logical charge (\Q\) should scale in proportion to y/V. In- 



III. TOPOLOGY AND THE HADRON 
SPECTRUM 

We address the question whether the statistics pre- 
sented in Fig. H suffices to expose a significant \Q\ depen- 
dency of hadronic states. To minimise statistical errors 
we subdivided each sample into two subsamples only, one 
containing configurations with \Q\ < 1.5 and one with 
\Q\ > 1.5. For the smaller volume at k = 0.1575 this 
division results in two subsamples of approximately 100 
configurations each. 

The most interesting test case is given by the effec- 
tive masses in the flavour singlet pseudoscalar channel, 
m°f (t), that should be particularly sensitive to the vac- 
uum topology. This is confirmed by the data displayed 
m Fig. |. We find the t)' mass on vacuum configurations 
with \Q\ < 1.5 to systematically lie below the correspond- 
ing mass obtained with the cut \Q\ > 1.5. Note that in 
the infinite volume limit one would expect such a sensi- 
tivity of correlation functions on \Q\ to disappear. For 
further details we refer to Ref. |E8| . 




t/a 

FIG. 4. The rj effective mass from a local-smeared corre- 
lation function on the 16 3 x 32 lattice at k — 0.1575, sepa- 
rately determined on configurations with topological charge 
\Q\ < 1.5 and \Q\ > 1.5. The horizontal error band is the 
fitted asymptotic mass obtained on the full sample p8[ . 

In Fig. ^ we show the corresponding flavour octet chan- 
nel quantities, i.e. n effective masses, where we anticipate 
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no such correlation with topology. This is indeed borne 
out by the data: not only at t — > oo but time slice by 
time slice no sensitivity to the value of the modulus of the 
topological charge is detected, with rather high statistical 
accuracy. The horizontal line with error band indicates 
the asymptotic large t results obtained from a fit to the 
entire data sample [E0|. 
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FIG. 5. The same as Fig. ^ for the 7r mass. 

Apart from the rf none of the standard mesonic, bary- 
onic and glueball-like states exhibit correlations between 
\Q\ and the respective effective masses. In Figs. ||- ||we 
illustrate this for the p, the nucleon and the A. In no 
case have we found any systematic effects on the static 
potential. Of course this does not exclude the possibil- 
ity of correlations between hadronic properties and the 
distribution of instantons and anti-instantons and their 
relative orientations in colour and position space [29] . 



IV. THE TOPOLOGICAL SUSCEPTIBILITY 

In Figs. |9|- [ll] we display three histograms of topologi- 
cal charge distributions. Each bin with width SQ, centred 
around Qj — jSQ, contains all measurements resulting in 
charges within the interval [Qj — 5Q/2,Qj + SQ / 2] , where 
SQ = 1 on the 16 3 x 32 lattice depicted in Fig. Hand 
SQ = 3 for the 24 3 x 40 lattices of Fig. [h] and Fig. y. In 
addition to the data we display Gaussian distributions, 



n(Q) 



S SQ 



: exp 



Q 2 



2(Q 2 



(10) 



In Fig. H we include the statistical uncertainties of the 
individual bins, An = 2Ti n t&' m cas\/'^; while in Fig. [l^ and 
Fig. [□] the error on the width of the distribution is re- 
flected by the error band around the central curve. 

The fitted topological susceptibilities \ = (Q 2 )/V as 
well as (Q) are displayed in Table [II. For orientation, 
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FIG. 6. The same as Fig. H for the p mass. 
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FIG. 7. Nucleon effective masses from smeared-smeared 
correlation functions. 
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FIG. 8. The same as Fig. for the A mass. 
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TABLE III. Topological susceptibilities. 



K 


L a 


(Q) 


X a 4 /10" 6 


X 1/4 r 


X 1/4 /MeV 


0.1560 


16 


0.04(32) 


69(6) 


0.466(11) 


183(4) (9) 


0.1565 


16 


0.48(48) 


83(13) 


0.505(22) 


199(8) (10) 


0.1570 


16 


0.60(34) 


57(10) 


0.476(22) 


187(8) (9) 


0.1575 


16 


0.05(33) 


41(7) 


0.472(21) 


186(8) (9) 


0.1575 


24 


0.84(1.33) 


44(7) 


0.479(20) 


189(8) (10) 


0.1580 


24 


1.92(1.01) 


31(12) 


0.465(44) 


183(17)(9) 


quenched 


16 


-0.68(33) 


80(14) 


0.505(22) 


199(9)(10) 




FIG. 9. Topological charge distribution at n — 0.156. 
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40 




FIG. 10. The same as Fig. H for k = 0.1575, L CT = 24 



we convert x into physical units in the last column of 
the table, using r^ 1 = (394 ±20) MeV. The distribution 
at k = 0.158 is not symmetric around zero anymore, as 
reflected by the value (Q) = 1.9 ± 1.0 and by Fig. [□]. 
Therefore, in this case, the resulting value of x should 
pass with a grain of doubt, notwithstanding the comfort- 
able error bars. A comparison between the L a = 16 and 
L a = 24 results at k = 0.1575 reveals that the level of 




FIG. 11. Topological charge distribution at k = 0.158. 

finite size effects on x is below the statistical errors. 

The Gell-Mann-Oakes-Renner (GMOR) relation con- 
nects the pion mass to the chiral quark mass, m, via the 
pion decay constant^ f v in the chiral limit and the chiral 
condensate, — S: 



j-2 2 



2mS + 0{m 2 



ill) 



/tt from the ap- 
yields V2f n = 



A direct lattice determination of 
propriate current matrix element 
Z A h\A 4 \Tr)/m^ -> 0.0496(34)a-\ in the limit ml -> 
[ p0[ . Allowing for an additional systematic error of 
10 % on to account for the fact that we have only de- 
termined the axial vector renormalization constants Za 
perturbatively, we arrive at the lattice estimate 

U = (0.238 ± 0.030) r 1 = (94 ± 13) MeV (12) 

at k = k c , using the value r f ^ 1 a = 0.1476^gQ [^TJ. In 
combining Eq. ( |ll| ) with Eq. (Q) one expectsP] 



X 



2n f 



0(mi 



(13) 



for small pseudoscalar masses and large Leutwyler- 
Smilga parameters, x = VmT, 3> 1 [ ^lp2] |. 

In Table IV we display the spatial lattice extents, 
Lrjarn-n as well as x and (Q 2 ) — \V- We estimate 
the (scheme independent) combination toE that ap- 
pears within x from Eqs. ( |IT|) - (|l2]). The topological 
charge fluctuations (Q 2 ) correspond to the average num- 
ber of non-interacting instantons that can be accomo- 
dated within the simulated lattice volume 33 3l] . It has 



2 Our normalisation conventions are such that / x ~ 92 MeV. 

3 This formula as well as the GMOR relation Eq. (0) only 
apply to nj > 2 since for iif = 1 no light pion exists (for a 
finite number of colours, N c < oo). 
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TABLE IV. The spatial lattice extents in units of the 
pion correlation length, the Leutwyler-Smilga parameters, 
x = VmT. and topological charge fluctuations (Q 2 ). 





V/a 4 


m^aha 


X 


(Q 2 ) 


0.1560 


16 3 


x 32 


7.14(4) 


16.0 ±4.1 


9.1 ±0.8 


0.1565 


16 3 


x 32 


6.39(6) 


12.9 ±3.2 


10.9 ± 1.7 


0.1570 


16 3 


x 32 


5.51(4) 


9.6 ±2.4 


7.5 ±1.3 


0.1575 


16 3 


x 32 


4.50(5) 


6.4 ± 1.6 


5.4 ±0.9 


0.1575 


24 3 


x 40 


6.65(6) 


26.1 ±6.6 


24.1 ±3.9 


0.1580 


24 3 


x 40 


4.77(7) 


13.4 ±3.4 


17.2 ±6.6 



0.08 




quenched 



FIG. 12. The topological susceptibility in physical units 
as a function of m£. The right-most data point is from the 
quenched reference simulation and has not been included into 
the fits. The dashed error band corresponds to the small mass 
expectation from the independently determined pion decay 
constant f„ = 0.238(30) r^ 1 . 



been argued j34j that x > 10 is already a large volume 
in the Leutwyler-Smilga 31 sense, such that Eq. Jl3| ) 
should safely apply to all but the 16 3 x 32 lattices at 
K = 0.157 and n = 0.1575. Even at the latter, more 
critical, K-value, where we have two lattice volumes that 
correspond toi~6 and x ~ 25 at our disposal, we do 
not find any volume dependence of \- 

In the regime of large one would expect \ m units 



of some reference scale like r$ to approach the quenched 
value as a smooth function of m n . Prior to the compar- 
ison with Eq. (H) one should in principle extrapolate 
the lattice results on \ to the continuum limit, for in- 
stance along a trajectory of constant m^ro. Having only 
= 5.6 data at our disposal we cannot yet perform this 
extrapolation but conjecture that our results are already 
sufficiently close to the continuum limit for Eq. (|l^) to 
hold. This assumption is plausible in view of the fact 
that the topological charges obtained from the gluonic 
and fermionic definitions agree with each other reason- 
ably well (with a renormalization constant Zp of order 
one) and of our previous studies of light hadronic quan- 



tities |2J|3QJ. 

For the purpose of the present investigation we shall 
consider the dimensionless topological susceptibility, x r o, 
as a function of the dimensionless combination to^Tq in 
order to carry out the chiral extrapolation. The results 
are plotted in Fig. |lj, together with the region that is 
allowed for by the lea ding order expectation, Eq. (|i~3]), 
for the /,r value of Eq. ([fj) . All unquenched data except 
for the point at the heaviest quark mass {m^r^ « 5) 
are consistent with this expectation. Note that the most 
chiral data point stems from the run with k = 0.158 
which, as can be seen from Figs. || and [fl], is n °t °f 
sufficient statistical quality for the purpose of topological 
studies. The right-most entry in Fig. [l^ is the result of 
our quenched reference study. 

A linear fit, excluding the largest mass point, to the 
parametrisation ( |l3|) renders — 0.255(ll)rQ 1 with 
X 2 /Ndf = 3.93/4. This value compares reasonably well 
with the expectation, Eq. (|l2[). On the other hand our 
data are consistent with a mass independence of x too: 
fitting them to a constant yields x r o — 0.0503(23) with 
X 2 /Ndf = 2.67/5, in agreement with the quenched ref- 
erence point, x r o = 0.065(11). 

Our present study is based on data obtained at a fixed 
value of (3 and the quark mass is varied by just tuning k. 
It is worthwhile to relate our findings to those recently 
presented by the UKQCD collaboration^] jl^] who simul- 
taneously vary (3 and k to change the quark mass while 
keeping the lattice spacing a fixed in units of ro . Apart 
from this they dispose of ensembles of gauge field con- 
figuration of comparable statistical sizes, using slightly 
coarser lattice spacings with an improved fermionic ac- 
tion. 

In order to facilitate this comparison we have super- 
imposed both data sets^] in Fig. |l3|. The data do not 
clearly contradict each other, however, we are unable to 
confirm the decrease of the topological susceptibility that 
has been reported in Ref. |15|, at pion masses within the 
range of present-day lattice simulations. 



V. CONCLUSIONS 

We have demonstrated that Wilson fermions are suit- 
able for lattice studies of topology related effects in the 
QCD vacuum. We have found agreement between the 
topological charge computed after cooling from the field 
strength tensor, Q = (F,*F)/(l6ir 2 ), and the fermionic 



4 Unfortunately, none of the other collaborations have con- 
verted their results into units of tq. 

5 The statistical errors quoted by the UKQCD collaboration 
are somewhat smaller than ours since no signs of autocorre- 
lation effects have been found in their study fjij. 
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FIG. 13. Our results in comparison to those reported by 
the UKQCD Collaboration The horizontal line is a con- 
stant fit to our data. The dashed error band corresponds to 
the small mass expectation from the pion decay constant. 



definition given by Smit and Vink with a renormaliza- 
tion constant Zp of order one. The trace has been esti- 
mated from diagonal improved Z^ noisy sources. As the 
next step we plan to investigate correlations between the 
topological charge density distribution and eigenvectors 
of 75M. 

The topological susceptibility at (3 = 5.6 has been 
found to be consistent with that of quenched studies at 
large sea quark masses. Our data at smaller mass val- 
ues is consistent with the asymptotic slope in that 
is expected from our independently determined n/ — 2 
pion decay constant. However, we are unable to unam- 
biguously confirm the decrease that has recently been re- 
ported by the UKQCD collaboration [jl5[ with a different 
fermionic action. Results obtained by the CP-PACS col- 
laboration |l4|] who cover a similar range of quark masses 
on somewhat coarser lattices do not clearly show this ten- 
dency either. We plan to clarify this issue in simulations 
at different (3 values which will enable us to perform a 
continuum limit extrapolation, and at smaller sea quark 
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